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(54) Adaptive equalization method and adaptive equalizer 



(57) During reception of a training signal in a re- 
ceived signal R(n) an estimated impulse response value 
H m (n) of an M-channel channel, and a tap coefficient G 
(n) of a linear filter 111 is calculated by an adaptive al- 
gorithm through the use of the received signal R(n) and 
the training signal b(n). For an information symbol of the 
received signal R(n) : the received signal R(n) is subject- 
ed to linear filtering with the most recently calculated tap 
coefficient G(n), and the linear filtering output 2(n) and 
the most recently estimated impulse response value H m 
(n) are used to calculate a soft decision value a. v In the 
second and subsequent rounds of equalization, the like- 
lihood b'(n) of a soft decision value A^[b(n)] from a de- 
coder is calculated, and a replica is generated by linear- 
filtering the likelihood b'(n) with an estimated impulse 
response value vector H L (n) obtained by approximating 
intersymbol interference with the current code b(n) to 
zero. A difference signal R c (n) between the replica and 
the received signal is calculated, and the estimated im- 
pulse response value vector H L (n) is used to update the 
tap coefficient G(n). Then the signal Z(n) is obtained by 
linear-filtering the difference signal R c (n) with the updat- 
ed tap coefficient G(n), and the signal Z(n) and the es- 
timated impulse response value vector H L (n) are used 
to calculate the soft decision value X. 1 . 
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Description 

BACKGROUND OF THE INVENTION 



5 [0001] The present invention relates to an adaptive equalization method and an adaptive equalizer for use : for ex- 
ample, in mobile communications to reduce the influence of fading. 

[0002] In a mobile communications environment signals transmitted from mobile stations arrive at a base station 
from various directions under the influence of surrounding buildings, trees and so forth. The influence of intersymbol 
interference, which results from the arrival of paths of large delay time differences in such a multipath channel, poses 

10 a serious problem for mobile communications. A known technique for lessening the influence of a signal waveform 
distortion by the intersymbol interference is an equalization technique that estimates the characteristic of the channel 
and performs equalization of delayed paths based on the estimated characteristic. One of equalization algorithms used 
in the equalization technique is a turbo-coding-technique-applied iterative equalization algorithm that achieves a very 
excellent error rate characteristic in low SNR. The turbo coding method implements, with practical coding-decoding 

15 complexity, an error rate characteristic that is almost a ceiling value obtainable with a channel coding method; only the 
presence of such an error rate is merely suggested in information theory. With the turbo coding method, the transmitting 
side transmits coded information sequence and its rearranged or interleaved version, and the receiving side repeatedly 
decodes them alternately with each other. Since the process of generating a coded sequence of convolution codes in 
the turbo coding is equivalent to the process of generation of the intersymbol interference in multipath channel, the 

20 decoding process of the turbo coding technique can be applied to the equalization of delayed paths. 

[0003] Fig. 1 illustrates in block form an example of a transmitter-receiver for the iterative equalization disclosed in 
Daryl Reynolds and Xiaodong Wang, "Low Complexity Turbo-Equalization for Diversity Channels," http://ee.tamu.edu/ 
reynolds/. A transmitter 10 encodes an information bit sequence d(i) by an encoder 11, then rearranges the coded 
output sequence by an interieaver 12, and modulates the carrier by a modulator 13, thereafter transmitting the mod- 

25 ulated output bit sequence b(i) over a channel 14. In a receiver 20 delayed paths are equalized by a SISO (Soft-Input 
Soft-Output) equalizer 21 for the received signal r(t), and the log-likelihood ratio (LLR) A^ between +1 and -1 for each 
code bit b(i) of the received signal r(t) is derived by the following equation (1). 



A,[bO)] = lo gpr[b( . ) = _ 1|r(t)] 



35 



p[r(t)|b(i) = +lj Pr[b(i) = +lj 
p[r(t)|b(i) = -l] *Pr[b(i)=-l] 

= A.,[b(i)] + MKi)] (1) 



[0004] In Eq. (1) the first term, A. 1 [b(i)] s is extrinsic or external information derived by subtracting a priori information 
40 AJ b (')]: by a subtractor 26, from the log-likelihood ratio A^bO)] provided from the SISO equalizer 21 . The second term, 
A^[b(i)] ; is the a priori information that is a rearranged or interleaved version, by an interieaver 25, of X^bfn)] obtained 
by subtracting ^(n)], by a subtractor 24, from a log-likelihood ratio A 2 [b(n)) provided from a SISO decoder 23 in the 
previous iteration. The initial value of X^[b(n)] is zero. The a priori information ^[b-0)] is used as an a priori LLR value 
of a coded bit string {b(i)}. The extrinsic information is restored by a de-interieaver 22 to X^bfJ)], which is input as a 
45 posteriori information to the SISO decoder 23. The SISO decoder 23 uses the a posteriori information to derive a a 
posteriori LLR value A 2 of the coded bit string {b(i» by the following equation (2). 



h 2 [b{n)]=X 2 Mn)]^[b(n)\ (2) 

50 

In the second and subsequent iterations, A..,[b(n)] is subtracted by the subtractor 24 from the a posteriori LLR value A 2 
derived by the SISO decoder 23 to obtain extrinsic information XJbln)], which is provided as a priori information ^[b 
(i)]via an interieaver 25 to the SISO equalizer 21 . Such iterations of equalization and decoding update the LLR value, 
providing a more probable LLR value. That is, the first log-likelihood ratio A 1 is calculated by the first adaptive equal- 
55 ization of the received digital signal, and the first log-likelihood ratio A^ is decoded to obtain the second log-likelihood 
ratio A 2 , which is used to perform the second adaptive equalization of the received digital signal to obtain the first log- 
likelihood ratio A 1; which is then decoded, and the decoded value is output. 
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[0005] In the iterative equalization the SISO equalizer and the SISO decoder use a MAP (Maximum A-posteriori 
Probability) algorithm, Log-MAP algorithm based on the MAP algorithm, or Max-Log-MAP algorithm that has reduced 
the computational complexity of equalization and decoding through use of approximate calculation. Still another algo- 
rithm is SOVA (Soft-Output Viterbi Algorithm) that provides bit reliability information while performing optimum decoding 

5 on a sequence-wise basis. In the case of using any of these algorithms for equalization by the equalizer, however, the 
complexity of the algorithm exponentially increases with respect to the maximum number of delayed symbols to be 
considered in the equalizer, resulting in very high computational complexity. For example, when BPSK is used as the 
modulation system, assuming that the storage capacity of one receiving channel is 10 symbols (the maximum number 
of delayed symbols to be considered in the equalizer being 10 symbols), the number of states becomes as large as 

w 2 10 =1024, exceeding practical computational complexity. The computational complexity could be reduced by using a 
linear equalizer for the iterative equalization. 

[0006] Fig. 2 illustrates a typical example of the linear equalizer. In the illustrated equalizer sampled received signals 
are applied to a cascade-connection of delay elements (of t-sec delay) 31 to obtain signals delayed by i-sec intervals 
in a sequential order, which are weighted by a variable weighting circuit 32 having set therein a coefficient by a coefficient 
15 adaptation circuit 34 and combined by an adder 33 to implement equalization. As the delay time i, the sampling period 
T orT/2 is usually chosen. Considering the equalizer as one filter, the linear equalizer achieves equalization by imple- 
menting a filter characteristic convenient for identifying the received signal. This filter basically implements an inverse 
characteristic of a channel distortion. 

[0007] When such a linear equalizer is used for the above-described iterative equalization, it is necessary to derive 
20 a soft decision value as the equalizer output. In the afore-mentioned literature by Daryl Reynolds and Xiaodong Wang, 
there is set forth a scheme in which, by deriving the soft decision value through the use of the linear equalizer, the 
iterative equalization can be achieved with less computational complexity than in the case of using the conventional 
algorithm such as MAP, SOVA or the like. In the literature a description is given of the application of the linear equalizer 
assuming the following discrete-time model as a channel model. A received signal sample r(n) is expressed as follows: 



25 



30 



35 



45 



50 



55 



J-l 

r(n) = £b(n - j)h(n; j) + ov(n), n=0, 1, N-l (3) 



v(n)=[v 0 (n)v 1 (n)...v M . 1 (n)] T (4) 



h(n;j) =[h 0 (n;j) h^nj) ...h M _., (n;j)] T (5) 



In the above, M is the number of channel outputs and 

r ( n )=[ r o( n ) r i( n )--- r M-i( n )] T is a vector representing a received sample value in each channel at time n, where T indicates 
40 a transposed matrix. For example, in the case of using a plurality of antennas, M corresponds to the number of antennas. 
B(n) represents a coded bit, and v(n) in Eq. (4) is a complex vector of an average 0, representing noise. In Eq. (5), h 
(n;j) is a vector representing the tap weighting coefficient of the channel, and J is the channel ISI length. Now, a received 
signal vector R(n), a noise vector V(n), a coded bit vector B(n) and a channel matrix H m (n) are defined as follows: 



R(n) = 



r(n + J-l) 
r(n) j 



(6) 

MJxl 
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V(n) = 



v(n + J-l) 
v(n) 



(7) 



MJxl 



10 



B(n)=[b(n+(J-1))b(n+(J-2))...b(n+1)b(n)b(n-1)...b(n-(J-1))] T 



H(n)=[h(n,0)...h(n;J-1)] 



(8) 
(9) 



15 



20 



H m (n) = 



H(n) 0 — 0 

0 H(n) \ 
\ 0 

0 ... 0 H(n) 



MJx(2J-l) 



25 



30 



35 



h(n;0) h(n;l) — h(n;J-l) 0 -•• 0 

0 h(n;0) h(n;I) •■- h(n;J-l) 0 0 

0 0 b(n;0) h(n;I) ... b(n;J-I) 0 

0 - 0 h(n;0) h(n;l) - h(n;J-l) 



00) 



JMJx(2J-I) 

[0008] The vectors R(n) and V(n) given by Eqs. (6) and (7) are column vectors of a matrix with MJ rows and one 
column, that is, MJ-dimensional column vectors. B(n) is a (2J-1)-dimensional column vector, H(n) is a matrix with M 
rows and J columns : and H m (n) is a matrix with MJ rows and (2J-1) columns. From the above, Eq. (3) becomes as 
follows: 



40 



R(n) =H m (n)B(n) +oV(n) 



(11) 



[0009] Next, the likelihood of the code bit string is defined by the following equation using the a priori LLR value 
[b(k)] of the code bit string b(i). 



50 



b'(k) =tanh[X 2 [b(k)]/2], n-(J-1) ^ k ^ n+(J-1) 



(12) 



[0010] Eq. (12) has the property that it rapidly approaches 1 or -1 upon deviation of a variable x of a hyperbolic 
function tanhx from 0, and the likelihood b'(k) represents an estimated value of b(k). Further, that estimated value 
vector of the code bit in the received signal vector R(n) which affects, as intersymbol interference, code bit b(n) at time 
n is given by the following equation. 



B' (n)=[b'(n+(J-1 ))b'(n+(J-2))...b'(n+1 )0b'(n-1 ). ..b' (n-(J-1 ))] T 



(13) 



In Eq. (13).. since B'(n) is an estimated value vector that is used to subtract from the received signal vector R(n) the 
intersymbol interference component for the code bit b(n) at time n, the element corresponding to b(n) of the vector B' 
(n) is set to zero. A difference vector R c (n) by subtracting an intersymbol interference component replica from the 
received signal vector R(n) is defined by the following equation using the estimated intersymbol interference code 



4 



3NSDOCID: <EP 1217798A2_L> 



1 



EP 1 217 798 A2 

vector B'(n). 

R c (n) = R(n)- H m (n)B'(n) = H Jn)[B(n)-B'(n)]+av<n) (14) 

[0011] The difference vector R c (n) is subjected to filtering expressed by the following equation through the use of a 
linear filter of tap coefficient G(n). 

Z(n) = G(n) H R c (n) (15 ) 

where H is a conjugate transposed matrix. The tap coefficient G(n) is so chosen as to minimize the mean square error 
between the code bit B(n) and the filter output Z(n) as expressed by the following equation. 

G(n)=arg min E^b^-G^R^n)!! 2 } 

=arg min G(n) H E{R c (n)R <: (n) H }G(n)-2 GCn^EfbCn)^)} (16) 

20 where arg min E{} represents that one of all J-dimensional vectors G whose content of {} is minimum. 
[0012] Expansion of first and second terms in Eq. (16) gives 
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E{R c (n)R c (n) H hH m (n)A(n)H m (n) H + o 2 l (17) 
E{b(n)R c (n)}=H m (n) ej (18) 

where I is a unit matrix and 

ACn^CovfBO^-B'Oi)} 

=diag[l-b ,2 (n+(J-l)),...,l-b' 2 (n+l),l,l-b ,2 (n-l), 

...,l-b' 2 (n-(J-l))] (19) 

In the above, Cov represents a covariance matrix and diag represents a diagonal matrix, ej is a vector of a length 2J- 
1 , a J-th element being 1 and other elements being 0s. By solving the normal equation with the gradient vector on the 
right-hand side of Eq. (16) to 0 and substituting Eqs. (17) and (18), the optimum tap coefficient G(n) is given as follows: 

G(n)=[H m (n)A(n)H m (n) H + a 2 ir 1 H m (n)e J (20) 
Substitution of the thus derived tap coefficient G(n) into Eq. (15) gives 

Z(n)=e J T H m (n) H tH m (n)A(n)H m (n) H +o 2 l]' 1 [R(n)-H m (n)B(n)] (21) 

In this case, since it is necessary that H m (n) be a full column rank, the antenna number M needs to be at least 2 or more. 
[0013] Assuming that the filter output can be approximated by a Gaussian distribution, 2(n) by Eq. (21) can be de- 
scribed as follows (see V. Poor and S. Verdu, "Probability of Error in MMSE Multiuser Detection," IEEE Trans. Infor- 
mation Theory., vol. IT-43, No. 3, pp. 858-871, May 1997). 

Z(n)=n(n)b(n)+n(n) (22) 
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where \i(n) is an equivalent amplitude of the output signal and i\{n) represents a Gaussian distribution of an average 
0 and the variance v 2 (n). Therefore, u.(n) and v(n) can be expressed as follows: 

ji(n)=E{Z(n)b(n)} 

=rH n ,(n) H [H m (n)A(n)H m (n) H W J r|- 1 H m (n)] w (23) 

10 where []j j indicates the element at the intersection of the J-th row and the J-th column of the matrix. 

v 2 =var {Z(n)}=n(n)-u 2 (n) (24) 
15 From the above, the extrinsic information derived by the linear equalizer can be deduced from the following equation. 

X rwrrtl - loo P(Z(n)lb(n)= + 1) ^ 4Real{Z(n)} 

^ [b(n)1 - 109 p(Z(n)lb(n)=-1)- 1-u<n) (25) 

20 [0014] To obtain the optimum tap coefficient G(n) with the above-described method, however, it is necessary that 
the inverse matrix calculation expressed by the following equation be conducted in Eq. (20) for each point in time; 
hence, much time is required for calculation. 

25 O(n) = [HJn)A(n)H m (n) H + oV (26) 

[0015] With the conventional method for iterative equalization by the linear equalizer, the inverse matrix calculation 
needs to be conducted for each point in time to update the tap coefficient-this gives rise to the problem of high com- 
putational complexity. 
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SUMMARY OF THE INVENTION 

[0016] It is therefore an object of the present invention to provide an adaptive equalization method and apparatus 
that permit reduction of the computational complexity for deriving the tap coefficient of the linear equalizer for iterative 
equalization and hence allow shortening the time for computation. 

[0017] According to a first aspect of the present invention, a received signal sample value and a known signal are 
used to estimate the impulse response of the channel, and a received signal sample value sequence is used to adap- 
tively calculate the tap coefficient of a linear filter. The tap coefficient is used to perform linear filtering of the received 
signal sample value sequence, and the result of linear filtering and the estimated value of the impulse response are 
used to derive a soft decision value. 

[0018] With this construction, in the case of using the linear equalizer for iterative equalization, the received signal 
sample value sequence is input thereto to adaptively derive the tap coefficient of the linear filter— this permits reduction 
of computational complexity as compared with the prior art method involving the inverse matrix calculation. 
[0019] According to a second aspect of the present invention, in second and subsequent processing by an adaptive 
iterative equalizer for the iteration of equalization and decoding, an estimated value vector of the code bit based on 
the likelihood of a soft decision value from a decoder is subjected to linear filtering with an estimated impulse response 
value H m (n) to generate a replica of the received signal. The replica is subtracted from the received signal to create a 
difference signal R c (n) having removed therefrom the intersymbol interference component, then the tap coefficient of 
the linear filter is calculated, for example, using the Matrix Inversion Lemma, from the estimated impulse response 
value vector H L (n) approximated with the intersymbol interference component regarded as zero, then the difference, 
then the difference signal R c (n) is subjected to linear filtering using the tap coefficient, and the result of linear filtering 
and the estimated impulse response value vector are used to derive the soft decision value. 

[0020] With this construction, the tap coefficient of the linear filter in second and subsequent iterative equalization 
by the linear equalizer can be calculated, for example, using the Matrix Inversion Lemma, by which it is possible to 
achieve adaptive equalization with less computational complexity than in the case of the conventional that involves 
the inverse matrix calculation of Eq. (20). 

[0021] According to a third aspect of the present invention, a received signal sample value sequence and a known 
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signal are used to estimate the impulse response of the channel, then a threshold value is set accordingly and com- 
ponents of values smaller than the threshold value are removed from the estimated impulse response value vector, 
that is, the dimension of the estimated impulse response value vector is reduced. Further, components corresponding 
to those of the values smaller than the threshold value and hence removed from the estimated impulse response value 
5 vector are removed from the tap coefficient of a linear filter calculated from the received signal sample value sequence 
and the received signal vector to thereby reduce their dimension as well. The soft decision value is derived using the 
estimated impulse response value vector, the tap coefficient of the linear filter and the received signal vector of the 
thus reduced dimensions. 

[0022] The above scheme permits reduction of the computational complexity for deriving the soft decision value 

10 when the linear equalizer is used for iterative equalization. 

[0023] According to a fourth aspect of the present invention, the likelihood of the soft decision value from a decoder 
is calculated in second and subsequent processing by an adaptive equalizer for iteration of equalization and decoding, 
then the tap coefficient of the linear filter is calculated, for instance, using the Matrix Inversion Lemma from the estimated 
impulse response value vector reduced in dimension in the first round of iteration and stored in a memory, then the 

is tap coefficient is used to perform linear filtering of a difference signal having subtracted therefrom an intersymbol 
interference component, and the result of linear filtering and the estimated impulse response value vector are used to 
derive the soft decision value. 

[0024] This scheme allows reduction of the computational complexity for deriving the soft decision value in the second 
and subsequent iterative equalization by the linear equalizer as well. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 



[0025] 



25 P'S 1 is a block diagram of a transmitter-receiver of the iterative equalization scheme; 

Fig. 2 is a block diagram of a linear adaptive equalizer; 

Fig. 3 is a block diagram illustrating the functional configuration for first equalization processing in a first embodi- 
ment of the adaptive equalizer according to the present invention; 

Fig. 4 is a block diagram of the functional configuration in second and subsequent rounds of equalization: 
30 Fig. 5 is a graph for explaining the effect of the present invention; 

Fig. 6 is a diagram depicting an example of a multipath channel environment; 

Fig. 7 is a block diagram illustrating the functional configuration for first equalization processing in a second em- 
bodiment of the adaptive equalizer according to the present invention; 

Fig. 8 is a block diagram of the functional configuration in second and subsequent rounds of equalization; 
55 F '9- 9 is a graph showing an example of an impulse response contained in the received signal; 

Fig. 10 is a graph showing an example of selecting effective paths from the impulse response; 

Fig. 11 is a graph for explaining the effect of the present invention; and 

Fig. 12 is a graph for explaining a channel model used for computer simulations. 

40 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 



EMBPDIMENT 1 



[0026] In the method set forth in the afore-mentioned literature that employs a linear equalizer for iterative equaliza- 
tion, the likelihood b'(k) of the code bit is calculated by Eq. (12) using the a priori LLR value >^[b(k)] of the code bit 
string. In the first equalization processing, since the a priori LLR value A^bfk)] in Eq. (12) is zero in the absence of the 
a priori information, the likelihood b'(k) is all zeros as shown below. 



so b'(k) = [00...0...] T 

Therefore, the difference vector R c (n) by Eq. (12) becomes as follows: 

55 R c ( n ) = R ( n ) -H m (n)B'(n) =R(n) (27) 



and A(n) becomes a (2J-1)x(2J-1) whose diagonal components are all ones as shown below. 
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A(n) = Cov{B(n)-B'(n)}=diag{1 ,1 ,...,1 , ...,1 ,1} 



(28) 



[0027] Accordingly, the optimum tap coefficient G(n) of the linear filter by Eq. (20) can be derived by the following 



From Eq. (29) the tap coefficient G(n) is dependent on the value H m (n) corresponding to the impulse response of the 
channel. In the case where temporal variations of the channel due to high transmission rate or the like, that is, where 
the duration of a received burst path is so short that the channel characteristic of the received signal does not change 
in that time interval, the tap coefficient G(n) is substantially constant. Accordingly, by precalculating the tap coefficient 
G(n) prior to the processing of the information symbol period and using the precalculated value for processing of each 
information symbol, the computational complexity can be reduced as compared with that in the conventional system 
that derives the tap coefficient G(n) for each information symbol. Moreover, since the above-mentioned temporal var- 
iations need not be taken into account as well in the calculation of the tap coefficient by Eq.(20), an adaptive algorithm 
can be used-this avoids the necessity for conducting the inverse matrix calculation to derive the tap coefficient G(n), 
permitting reduction of computational complexity. 

[0028] In view of the above, in the first embodiment of the present invention, the tap coefficient calculation is con- 
ducted by an adaptive algorithm to obtain G(n) in a known signal period or training signal period prior to the processing 
of the information symbol period, and the value thus obtained is used to perform the processing in the information 
symbol period. This eliminates the needs for the inverse matrix calculation and the derivation of G(n) at each point in 
time of the information symbol, enabling the computational complexity to be made lower than in the conventional 
scheme. 

[0029] Figs. 3 and 4 depict examples of functional configurations for first processing and the second and subsequent 
processing for equalization by the adaptive equalizer according to the first embodiment of the present invention. The 
illustrated adaptive equalizer corresponds to the equalizer 21 in Fig. 1 . Fig. 3 shows the functional configuration for 
the first equalization processing, and Fig. 4 for the second and subsequent equalization processing. The broken line 
in Fig. 4 indicates the portion that is not used at the current point in time (in second and subsequent equalization 
processing) but used in the first equalization processing. 

[0030] The adaptive equalizer of the first embodiment comprises an impulse response estimating part 100, a linear 
filter part 110, a soft decision value calculating part 120, a switch 101, a likelihood calculating part 130, and a switch 
160. The impulse response estimating part 100 is made up of a transversal filter 103, a subtractor 1 05 : and an adaptive 
estimating part 106. The linear filter part 110 is composed of a transversal filter 111 , a tap coefficient calculating part 
112, and a subtractor 113. 

[0031] A description will be given, with reference Fig. 3, of the first equalization processing in the iterative equalization. 
In this instance, the switch 1 01 is connected to a known signal, for example, training signal input terminal 1 02, and the 
switch 160 is connected to an input terminal 104. In the impulse response estimating part 1 00, a training signal sequence 
b(n), provided as a known pattern from the training signal input terminal 102, and a difference value (vector) R c (n), 
obtained by subtracting an output replica H m (n)b(n) of the transversal filter 103 by the subtractor 105 from a received 
signal sample value sequence R(n) from the input terminal 104 are used to estimate the impulse response H m (n) of 
the channel in a calculating part 1 06A of the adaptive estimating part 1 06. This impulse response estimating calculation 
needs only to minimize the square of R c (n) through the use of an adaptive algorithm utilized in an ordinary adaptive 
equalizer. The estimated impulse response value H m (n) having converged at the point of termination of the training 
signal sequence b(n) is stored in a storage part 106B. 

[0032] In the linear filter part 110, the received signal sample value sequence R(n) is input to the transversal filter 
111 as a linear filter, by which it is subjected to linear filtering with the tap coefficient G(n) calculated in the tap coefficient 
calculating part 112, then the filtered output Z(n)=G H (n) H R(n) is subtracted by the subtractor 113 from the training 
signal sequence b(n), and the difference is input to the tap coefficient calculating part 112 to update the tap coefficient 
G(n) of the linear filter. This updating calculation utilizes an adaptive algorithm. That is, in the first iteration of equali- 
zation, since R(n)=Rc(n) from Eq. (27), the tap coefficient G(n) is so determined by the adaptive algorithm as to minimize 
||b(n)-G(n) H R(n)|| from Eq. 16. On this account, the tap coefficient G(n) can be calculated with less complexity than in 
the case of conducting the inverse matrix calculation by Eq. (20). The updating of the tap coefficient G(n) stops with 
the termination of the training signal sequence b(n), and the tap coefficient at this time is held unchanged. 
[0033] The received signal after termination of the training signal sequence b(n), that is, the information-symbol 
received signal sample sequence R(n), is subjected to linear filtering by the transversal filter 111 with the most recently 



5 



calculation. 



G(n)= [H m (n)H m (n) H + a 2 ir 1 H m (n)e J 



(29) 
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updated tap coefficient G(n). The filtering output 2(n) and the ultimately estimated channel impulse response value 
H m (n) are used to calculate, as the first equalization output, the soft decision value A^bfn)] in the soft decision value 
calculating part 120. In the first iteration of equalization, since the a priori information X^n^O, A 1 =fb(n)]=A. 1 [b(n)]. 
The soft decision value ^[b(n)] is fed to the de-interleaver 22 (see Fig. 1). 

5 [0034] Next, a description will be given first of the principle of equalization in the second and subsequent rounds of 
iterative equalization and then of concrete equalization processing with reference to Fig. 4. In the second and subse- 
quent rounds of equalization processing, the received signal R(n) used in the first iteration of processing is input again. 
On the other hand, when the likelihood b»(n) is calculated by Eq. (12) from the output of the decoder 23 (Fig. 1 ) produced 
for the result of the first iteration equalization, to be more precise, from the soft decision value ^^(n)] of the code bit 

io fed from the interleaver 25, the likelihood b'(n) is approximately +1 or -1 if Eb/No (Bit energy-noise ratio) is high. In 
other words, the estimated intersymbol interference code vector of Eq. (1 3) can be approximated as expressed by the 
following equation. 

15 B'(n)-[±1 : ±1, ±1,0, ±1,..., ±1, ±1] T (30) 

The use of thus approximated b'(n) permits approximation of Eq. (19) to the following equation. 
20 A(n)=Cov{B(n)-B'(n)}:=diag{0,0 0,1 ,0,... ,0,0} (31) 

Substituting Eq. (31 ) into Eq. (20), the tap coefficient G(n) of the linear filter can be obtained by conducting the following 
calculation. 



25 



30 



G(n)=[H L (n)H L H (n) + a 2 l]' 1 H L (3 2) 



where H L (n) is the vector of one column of M by J rows expressed by 

H L (nMh(n; J-1)...h(n; 0)] T =[h 0 (n; J-I^.h^n; J-1) 
h 0 (n; J-^.h^n; J-2)...h 0 (n; 0)...h M ^(n; 0)] T 



[0035] In this instance, the use of the Matrix Inversion Lemma in the calculation of the tap coefficient G(n) of the 
linear filter by Eq. (32) allows reduction of the computational complexity therefor. With such approximation of the es- 
timated intersymbol interference code vector B(n) as expressed by Eq. (30), H L (n) becomes the vector of one column 
of M by J rows as mentioned above -this avoids the need for conducting the inverse matrix calculation, permitting 
substantial reduction of the computational complexity as compared with that needed by the conventional scheme: 
Thus, according to the present invention, in the second and subsequent rounds of iterative equalization the computa- 
tional complexity in the tap coefficient calculating part for deriving the tap coefficient G(n) of the linear filter can be 
reduced through utilization of the approximation by Eq. (31) and the Matrix Inversion Lemma. 

[0036] Fig. 4 illustrates in block form an example of the functional configuration in the second and subsequent rounds 
of equalization by the adaptive equalizer of the first embodiment. The change-over switch 101 is connected to the 
output of the likelihood calculating part 130, and the switch 160 is connected to the output of the subtractor 1 05. The 
soft decision value ^[bfn)] of each bit of the code sequence, which is the output from the decoder 23 (the output from 
the interleaver 25) in the previous round of iterative equation, is input to the likelihood calculating part 130 from a 
terminal 130, which conducts the calculation of Eq. (12), outputting the likelihood b'(n). The likelihood b'(n) is filter- 
processed by the transversal filter 103 using the estimated channel impulse response value H m (n) stored in the storage 
part 1 06B, obtaining H m (n)B'(n) in Eq. (1 4) as a replica of the intersymbol interference in the received signal R(n). The 
output replica of the transversal filter 1 03 is subtracted by the subtractor 1 05 from the received signal value sequence 
R(n) to derive the difference value R c (n) of Eq. (14). 

[0037] The tap coefficient G(n) of the linear filter is calculated by Eq. (32) using, as the input therefor, only the esti- 
mated channel channel impulse response value H L (n) stored in the storage part 1 06B unlike in the scheme for deriving 
the tap coefficient of the linear filter in the first iteration of iteration. 

[0038] This calculation utilizes the afore-mentioned Matrix Inversion Lemma. That is, when this Lemma is expressed 
by A=B' 1 +CD- 1 C H , where A and B represent (M by M) positive definite matrixes, C an (M by N) matrix, D an (N by N) 
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positive definite matrix, an inverse matrix of A is given by 



A* 1 =B-BC(D+C H BBC)" 1 C H B 



(33) 



5 



Applying this Lemma to the calculation of the inverse matrix in Eq. (32), 



H L (n)H L H (n)+o 2 l=B' 1 +CD" 1 C H ! 



10 



H L (n)H L H (n) =CD~ 1 C H , 



15 




h(n) =C, 



20 



H L H (n)=C ! 



H 



25 



By using them to calculate Eq. (33), the inverse matrix calculation in Eq. (32) is solved. Incidentally, Eq. (33) also 
contains an inverse matrix calculation (D+C H BBC)* 1 , but this inverse matrix is also scalar, and hence it can similarly 
be calculated. 

[0039] The thus calculated tap coefficient value G(n) is used to perform linear filtering of the difference value R c (n) 
by the transversal filter 111, and the filtering result Z(n) and the estimated channel impulse response value H L (n) stored 
in the storage part 106B are used to calculate Eq. (25) in the soft decision value calculating part 120 to derive the soft 
decision value A 1 [b(n)]. The equivalent amplitude jx(n) in Eq. (25) is given by the following equation by substituting Eq. 
(31)forA(n) in Fig. (23). 



Incidentally, in the first iteration of equalization the tap coefficient calculating part 112 of the linear filter part 110 may 
also use the Matrix Inversion Lemma to calculate the tap coefficient G(n). 

[0040] As described above, this embodiment permits reduction of the computational complexity for deriving the tap 
coefficient of the linear equalizer that is used for iterative equalization. Fig. 5 shows an error rate characteristic obtained 
by computer simulations, for example, when an RLS (Recursive Least Squares) algorithm of excellent convergence 
characteristic is used as an adaptive algorithm for deriving the tap coefficient of a linear equalizer used for iterative 
equalization. 

[0041] In Fig. 5, the ordinate represents BER (Bit Error Rate) and the abscissa Eb/No (bit energy-noise ratio). The 
modulation system used was BPSK (two-phase shift), the information transmission rate was 1 2 Mbps, and the maximum 
Doppler frequency indicating the fading variation rate was 1000 Hz. The frame configuration used is a unique word of 
a known pattern for first 1 28 symbols, then the channel impulse response H m (n) was estimated, and 1 28 symbols were 
added as information symbols to the end of the training symbols. The RLS algorithm was used for the estimation of 
the impulse response. The channel was set as a five-wave Rayleigh fading channel of the same level, and the receiving 
side was to perform two-branch diversity and iterative equalization by a five-tap linear filter. Further, SOVA (Soft -Output 
Viterbi Algorithm) was used for the decoding in the iterative equalization, and the number of iterations was two. 
[0042] In Fig. 5, the curves AO and A1 show tap coefficients of the linear filter obtained by conducting the inverse 
matrix calculation at respective point in time in the cases where the number of iterations were one and two, respectively. 



35 



Kn) = [H L (n) H [H L (n)H L (n) H +o 2 I]- , H L (n)3 
=H L (n) H G(n) (34) 



40 
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The curves BO and B1 show the tap coefficients G(n) derived by the scheme according to the first embodiment of the 
present invention, the curve BO indicating the tap coefficient values by using the adaptive algorithm and the curve B1 
the tap coefficient values by using the Matrix Inversion Lemma when the number of iterations was two. In the case of 
the present invention, degradation of the bit error ratio is only on the order of 0.1 dB as compared with the ratio in the 
case of conducting the inverse matrix calculation at each point, it can be regarded as about the same as in the latter 
case. The present invention does not need to conduct the inverse matrix calculation expressed by Eq. (26), and hence 
it permits substantial reduction of the computational complexity for deriving the tap coefficient of the linear equalizer 
that is used for iterative equalization. Incidentally, it was confirmed that the same goes for the cases where the known 
pattern is 128-symbol and the number of information symbols is 256 and where the known pattern is 128-symbol and 
the number of information symbols is 512. This means that when the duration of a burst path (the frame length) is short 
about the same as these information symbol durations, impulse response variations of the channel can be ignored, 
and the bit error rate is approximately equal to that in the case of conducting the inverse matrix calculation at each 
symbol time, and with less computational besides. Letting M and J represent the number of antennas (channels) and 
the length of intersymbol interference (the number of taps of the linear filter 111), respectively, the computational com- 
plexity is around 0((MJ) 3 ) in the case of conducting the inverse matrix calculation at each symbol time, but in the above- 
described embodiment of the present invention the computational complexity is approximately 0((MJ) 3 ), which is sub- 
stantially lower than in the above case. 

[0043] While the examples of Figs. 3 and 4 are shown to include two transversal filters 103 and 111 with a view to 
facilitating a better understanding of the function configuration, the number of filters need not always be two but may 
also be one. Moreover, the tap coefficient calculating part 112 and the adaptive estimation part 106 are shown as 
different blocks, but they can be combined into one unit. 

[0044] As described previously, the one mode of the first embodiment utilizes the fact that the impulse response of 
the channel does not so vary over a short range of time, and hence it is characterized in that the training signal at the 
beginning of the received burst path (frame) or the known signal is used to obtain the estimated impulse response 
value H m (n), which is used to the subsequent processing. 

[0045] Further according to the other mode of this embodiment, in the second iterations of iteration, only a J-th 
element of a diagonal matrix is set to 1 which is a covariance matrix of the difference between the code bit vector B 
(n) and the estimated intersymbol interference code vector B'(n) as expressed by Eq. (31), and the other diagonal 
elements are approximated to zeros, by which the calculation for the tap coefficient G(n) becomes a vector calculation 
as indicated by Eq. (32)-this avoids the necessity of conducting the inverse matrix calculation and hence allows re- 
duction of computational complexity. In the first processing the known signal (training signal) need not always be used, 
but instead the output from a discriminating part 150, which discriminates the output from the transversal filter 103 
between two value in binary, may be used in place of the signal from the output of the switch 101 as indicated by the 
broken lines in Fig. 3. The estimated impulse response value H m (n) ultimately obtained by the first processing is used 
in the second and subsequent processing. 

EMBODIMENT 2 

[0046] In the first embodiment described above with reference to Figs. 3 and 4, all those of path components so- 
called multipath received signal R(n) arriving at the input terminal 104 via different paths which are within the range of 
the maximum number J of delayed symbols to be considered are equalized as will be seen from Eqs. (5), (6), (8) to 
(11), (15) and (25). According to the channel environment of the channel, however, path components of relatively large 
received power may sometimes arrive apart in time from other path components as indicated by multipath components 
d1, d2, ... in Fig. 6. In such a channel environment, too, the first embodiment performs equalization of all the received 
path components concerned. 

[0047] That is, in the first embodiment, path components within the range of the maximum number J of delayed 
symbols, even if very small, in extreme cases, zero in received power, are also subjected to equalization, and conse- 
quently, much time is needed for calculation. As a solution to this problem, this embodiment is adapted to omit particular 
ones of the dimensions corresponding to multipath components of negligibly small received power to reduce the sizes 
of vectors and matrixes to be handled, thereby reducing the computational complexity and hence shortening the cal- 
culating time. 

[0048] Calculating the likelihood b*(n) of the soft decision value A^bfn)] for the code bit string b(n) that is the decoder 
output for the result A 1 [b(n)] of the first iteration of equalization, by Eq. (12) as in the second iterations of equalization 
in the first embodiment, the calculated value is virtually +1 or -1 . Through utilization of this fact for the approximation 
of the estimated intersymbol interference code vector B'(n) as expressed by Eq. (30), the tap coefficient G(n) of the 
linear filter 111 can be calculated using the Matrix Inversion Lemma as follows: 
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G(n)= [H L (n)H L (n) H + a^'H^n) 

= H L (n)/(a 2 + H L (n) H H L (n)) (35) 

5 

where H L (n) is a column vector of a M by J length as shown below. 



H L (n) = [h 0 (n; J-1 )...□„., (n; J-1) h 0 (n; J-^.h^n; J-2) ... 



w 



.h 0 (n;0)...h M . 1 (n:0)] (36) 



By the approximation ot A(n) in Eq. (23) as indicated by Eq. (32), the equivalent amplitude p.(n) of the output signal is 
also calculated as expressed by (34). 
15 [0049] On the other hand ; Z(n) is given by 

Z(n)=G(n) H R c (n) (37) 

20 as is the case with Eq. (15). By using in Eq. (25) the linear filter output Z(n) obtained by Eq. (37), the extrinsic information 
value A-,[b(n)] derived by the linear equalizer can be used as the results of second and subsequent equalizations. 
[0050] Now, assume that path components j= 1 to j=J-2 within the range of the maximum number J of delayed symbols 
to be considered in the equalizer have not arrived but that path components j=0 and j=J-1 have arrived. If the impulse 
response of the channel can accurately be estimated, other path components than effective ones are such that h(n; 

25 1) : h(n;2) ; ...,h(n;J-2)=0. Accordingly, it follows from Eq. (35) that 

H L (nHh 0 (n;J-1)...h M . 1 (n;J-1)0...0 0...0 h 0 (n;0) 

& ...h^niO)} 7 (38) 

[0051] Since the tap coefficient G(n) of the linear filter can be calculated from H L (n) as shown by Eq. (35), it follows that 
35 G(n)= [g 0 (n;J-1 )...g M . 1 (n;J-1) 0...0 g 0 (n;0)...g M . 1 (n;0)] T (39) 

[0052] The derivation of the soft decision value A.,[b(n)] by the linear equalizer by Eq. (25) requires the calculation 
of the equivalent amplitude \i{n) of the output signal by Eq. (36) and the linear filter output Z(n) by Eq. (37). The 
equivalent amplitude p.(n) and the linear filter output Z(n) are calculated as shown below. 

40 

H(n) = H L (n) H G(n) 

= ho(n; J-l)xgo(n;J-l)+...+ h M .!(n; J-l)xg M _,(n;J-l)+...+0x0+_.. 

45 

. . .+0x0+. - •+ h o (n;0)x g0 (n;0)+. . .+ h M .,(n;0)xg M .,(n;0) 

(40) 

50 [0053] Here, in the case where zero-elements of H L (n) H and G(n) are omitted, that is, where the dimensions of H L 
(n) H and G(n) are lowered as shown below, too, the value of u.(n) naturally remains unchanged. 



55 



H L '(n) H = [h 0 (n;J-1)...h M . 1 (n:J-1) h 0 (n;0)...h M . 1 (n;0)] (41) 
G , (n) = [g 0 (n;J-1)...g M . 1 (n;J-1)g 0 (n;0)...g M . 1 (n;0)] T (42) 
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H(n) = H L Xn) H G'(n) 

= h 0 (n;J-l)x g0 (n;J-l>(-...+ h M -i(n;J-l)xg M _,(n;J-l) 

+ho(n;0)xgo(n;0)+. . .+h M -t(n;0)xg M .,(n;0) (43) 

And, Z(n) becomes as follows: 

Z(n) = G(n) H R c (n) 

= g 0 (n; J-l)x ro (n+J-l>J-...+ gM.,(n; J-l)xr M _,(n+J-l)+... 

+0xr 0 (n+J-2)+. . .+0xr M _,(n+J-2)+ +Oxr 0 (n+l )+. . .-H)x 

r M .,(n+l)+go(n;0)xro(n)+. ..+g M .,(n;0)xr M .,(n) (44) 

[0054J In Eq. (44), since G(n) H has O-elements, the values of R c (n) corresponding to the O-elements are unnecessary 
in the derivation of Z(n) as shown below. 

R c '(n) = [r 0 (n+J-1)...r M . 1 (n+J-1)r 0 (n)...r M . 1 (n)] (45) 

Z'(n) = G'(n) H IV(n) 

= go(n;J-l)x ro (n+J-l )+...+ g M .,(n; J-l)xr M _,(n+J-l)+... 
. . .+ g() (n;0)xr 0 (n)+. . .+ g M .,(n;0)xr M . 1 (n) (46) 

[0055] Since the dimension of R c (n) can be lowered as described above, the dimension of the received signal vector 
R(n) can also be lowered by Eq. (14). Further, since the dimension of the intersymbol interference component replica 
H m( n ) B '(n) can also be reduced, the dimension of the channel matrix H m (n) can be reduced as well. 
[0056] From the above, it can be seen that calculations for the path components having arrived can be dispensed 
with in the calculation for deriving the extrinsic information by the linear equalizer. That is, in the estimated impulse 
response value of the channel., the tap coefficient G(n) of the linear filter the received signal vector R(n) and R c (n) 
obtained by subtracting the intersymbol interference component from the received signal vector, those elements of 
vectors and matrixes corresponding to path components having not arrived can be removed. Accordingly, it is possible 
to decrease the number of dimensions of vector components and the number of elements of matrixes necessary for 
the calculation involved in the derivation of the soft decision value of the linear equalizer. 

[0057] However, the estimated impulse response value does not become zero due to an error in the estimation of 
the impulse response of the channel. Hence, a certain threshold value h lh is set for comparison and discrimination of 
the estimated impulse response value as described below. For example, if h(n;1 ),h(n;2),...h(n;J-2)<h th! then it is re- 
garded as 

h(n;1) J h(n;2),...h(n;J-2)=h th . Letting L represent the number of effective received path components larger than the 
threshold value h th , the dimension of the estimated impulse response vector H L (n), which indicates ail estimated impulse 
response values within the range of the maximum number J of delayed symbols, can be lowered from an M*J-th 
dimensional vector to an IvTL-th dimensional vector H' L (n) as expressed below. 

H L (n) =[h 0 (n; J-1)...h M . 1 (n; J-1) 0... 0 h o (n;0)...h M _ 1 (n;0)] T (47) 
H L '(n)=[h 0 (n; J-IJ-.h^ (n; J-1)h 0 (n;0)...h M _ 1 (n;0)] T (48) 



13 

BNSDOCID: <EP 121 7798 A2J_> 



10 



EP1 217 798 A2 

[0058] Further, the number of row elements of the channel matrix H m (n) can be lowered from 2J-1 to 2L-1 . 



H m (n) = 



h(n;0) h(n;l) •• h(n;J-l) 0 
0 h(n;0) h(n;l) — b(n;J-l) 



0 
0 



h(n;0) h(n;l) 
0 h(n;0) 



h(n;l) 



h(n;J - 1) 



0 
0 



b(n;J - 1) 



JMJx(2J-l) 



(49) 



15 



tt Y x [bOnO) 0 - h(n;J-l) 0 0 ] 

W 0 ... 0 h(n;0) 0 ... h^j.^J 



MLx(2J-l) 



(50) 



20 



25 then 



h(n;1),h(n;2),...h(n;J-2)<h th! 



9 0 ( n ; 1 ) ■ g M -i( n > 1 )'go( n i 2 ) -g M . 1 (n;2) ! ...,g 0 (n;J-2)...g M . 1 (n;J-2)- 0 

30 is set; the tap coefficient vector G(n) of the linear filter in the second iterations of iteration is given by the following 
equation. 



G'(n)= tgo(n;J-1)■•g M - 1 (n;J-1)go(n;0)...g M . 1 (n;0)] , 



(51) 



35 



40 



This is apparent from Eq. (35) in the calculation of H L , (n) H H L , (n). Moreover the tap coefficient vector G(n) can also be 
lowered from the M*J4h dimensional vector to the M*L-th dimensional vector G'(n). 

[0059] Similarly, the received signal vector R(n) and the difference vector R c (n) can also be lowered from the M*J- 
th to the M*L-th dimension. 

[0060] The equivalent amplitude u.(n) of the output signal can be calculated by Eq. (43) from H L '(n) and G'(n), and 
the linear filter output Z'(n) can be calculated by Eq. (46) from G'(n) H and R c '(n). Accordingly they can each be derived 
by calculations of vectors of low dimensionality. 

[0061] Furthermore, it is possible to lower the dimension of the number of column elements at positions of path 
components regarded as zeros in the channel matrix Hm'(n) of Eq. (50) generated from the estimated impulse response 
value as expressed by the following equation, and the dimension of the channel matrix for deriving the soft decision 
value by Eq. (25) can also be reduced. In an example in which only paths of j-0 and j=J are set as effective paths, the 
number of column elements can be decreased from (2J-1) to (2L-1). 



50 



H Yn> = r b(n;0) h(n;, " 1) ° 1 (52\ 
mKJ I 0 b(n,0) MKJ-ljJ^n 



MLx(2L-l) 

55 [0062] Figs. 7 and 8 illustrate in block form examples of functional configurations in the first iteration and second and 
subsequent rounds of equalization. In this embodiment, an identifying part 140 and a switch 170 are additionally pro- 
vided in the functional configuration of the first embodiment shown in Figs. 3 and 4, and the tap coefficient calculating 
part 112 is composed of a calculating part 11 2A, a storage part 112B and an identifying part 112C. The adaptive 
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20 



25 



30 



35 



40 



50 



estimation part 106 is also added with an identifying part 106C. 

[0063] In the first place, the switch 101 is connected to the terminal 102 : and the switches 160 and 170 are both 
connected to the terminal 104. In the first iteration of equalization the processes for obtaining the estimated impulse 
response value H m (n) and the tap coefficient G(n) are the same as in the first embodiment. Further, in Fig. 7, effective 
paths are identified by the identifying part 106C based on the estimated impulse response value H m (n) having con- 
verged at the end of the training signal. Assume, for example, that estimated impulse response values are obtained 
at respective symbol delay times within the range of the maximum number J of delayed symbols as shown in Fig. 9. 
The effective path identifying method is, for instance, to select paths above a certain threshold value h th as depicted 
in Fig. 10. The threshold value h th is for instance, the average power of all received paths, and received paths above 
the average power are selected. The effective paths thus identified are each represented, for example, by the number 
of delayed symbols and stored in the storage part 106B together with the corresponding tap coefficient value h. This 
data may be stored as data of the vector H L '(N) given by Eq. (48) or as data of the matrix H m *(n) given by Eq. (48); 
and it can easily be obtained from the one to the other and vice versa. Further, only those of the tap coefficient vectors 
G(n) of the linear filter 111 which correspond to the effective paths identified by the identifying part 106C are each 
stored as the tap coefficient vector G'(n) in the storage part 112B. 

[0064] Moreover, the switch 160 is switched to the identifying part 140, in which only an information symbol received 
signal R*(n) corresponding to the effective path in the signals received after the training signal, that is, in the information 
symbol received signals R(n), is selected, and the thus selected signal is subjected to linear filtering by the transversal 
filter 111 using its tap coefficient G'(n) stored in the storage part 112B most recently. The linear filtering output Z*(n) 
and the most recently estimated impulse response value H' L (n) are used to calculate the soft decision value A t [b(n)] 
in the soft decision value calculating part 120. 

[0065] In the second and subsequent rounds of equalization, the switch 1 01 is connected to the likelihood calculating 
part 130, the switch 1 60 to the identifying part 140 and the switch 170 to the output of the subtractor 105 as depicted 
in Fig. 8. The likelihood b'(n) representing the estimated intersymbol interference code, provided from the likelihood 
calculating part 130, is subjected to filtering by the transversal filter 103 using the estimated impulse response value 
H L (n) stored in the storage part 106B, and the output replica H L (n)b'(n) is subtracted by the subtractor 105 from the 
received signal sample value sequence R'(n) to obtain the difference value R c '(n). For only the path identified in the 
identifying part 106C of the impulse response estimating part 100, the difference value R c '(n) is identified by the iden- 
tifying part 140, and the identification result is input to the linear filter 111 . The tap coefficient G'(n) of the linear filter 
111 is derived from only the estimated impulse response value H L '(n) stored in the storage part 106B, by Eq. (35) 
through utilization of the Matrix Inversion Lemma, and the tap coefficient G'(n) is stored in the storage part 112B. The 
thus stored estimated impulse response value H L '(n), the tap coefficient G'(n) of the linear filter 11 1 and the identified 
difference value R c '(n) are used to calculate the soft decision value A 1 [b(n)] in the soft decision value calculating part 

[0066] As described above, the second embodiment permits reduction of the computational complexity for deriving 
the soft decision value of the linear equalizer. Fig. 11 shows the results of computer simulations of the bit error rate 
characteristic in the case of using the RLS (Recursive Least Squares) algorithm of excellent convergence characteristic 
for deriving the tap coefficient of the linear equalizer in the first iteration of iterative equalization and the Matrix Inversion 
Lemma in the second and subsequent iterations of iterative equalization. 

[0067] The modulation system used is BPSK; the information transmission rate is 10 Mbps; the maximum Doppler 
frequency indicating the fading variation is 1 000 Hz; and the frame configuration is one that a training signal of a known 
pattern of 63 symbols is added as an information signal to 256 symbols. A channel, in which only path components 
j=0 and j=4 as indicated by the solid lines in Fig. 12 had arrived, was assumed as a channel model. The receiving side 
performs two-branch diversity and uses SOVA (Soft-Output Viterbi Algorithm) for decoding. On the assumption that 
the channel was ideally estimated, the number of iterations was set to three. 

[0068] In Fig. 1 1 , the curves A1 , A2 and A3 show the results obtained in the cases where: equalization was performed 
once and iterated twice and three times, respectively; the tap number of the linear equalizer was considered for 0 to 
4T delayed paths; and the entire information of five taps was used as the tap coefficient of the linear equalizer. The 
curves B1, B2 and B3 show the results obtained in the cases where: equalization was similarly performed once and 
iterated twice and three times, respectively; and two-tap (only 0 and 4T delayed paths) information was used as the 
tap number of the linear equalizer. 

[0069] As is evident from Fig. 1 1 , a comparison of the results obtained with the use of five-tap and two-tap information 
as the tap number of the linear equalizer indicates that two and three iterations of equalization provide about the same 
error rate characteristic with low computational complexity. 

EFFECT OF THE INVENTION 

[0070] As described above, the first embodiment of the present invention enables the linear equalizer for iterative 
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equalization enables to reduce the computational complexity for calculating the tap coefficient of the linear filter and 
hence shorten the time for the calculation. 

[0071] According to the second embodiment of the invention, in a channel environment in which path components 
of relatively large power have arrived apart in time from other path components, it is possible to adaptively reduce the 
5 computational complexity for calculating the extrinsic information by the linear equalizer that performs iterative equal- 
ization. 



Claims 

10 

1. A method for adaptive equalization method of a received signal, said method comprising the steps of: 

(a) calculating an estimated impulse response value of a channel based on said received signal and a known 
signal, and storing said estimated impulse response value in a memory; 
15 (b) calculating the tap coefficient for linear filtering based on said received signal and said known signal; 

(c) linear-filtering said received signal with said tap coefficient; and 

(d) calculating a soft decision value from the result of said linearfiltering of an information signal of said received 
signal following at least said known signal and said stored estimated impulse response value. 

20 2. In decoding process for decoding a received code bit b(n) by iterating adaptive equalization of said received code 
bit b(n) a plurality of times, an adaptive equalization method in which a second soft decision value X^[b(n)] generated 
by the soft decision of said code bit b(n) is provided as a priori information, and the equalization of a received 
signal R(n) of an M-channel channel through utilization of said a priori information and outputting of a first soft 
decision value A^ttyn)] are iterated to adaptively equalize said received signal R(n) : said method comprising the 

25 steps of : 

in a first round of equalization, 

(a) calculating an estimated impulse response value H m (n) of said M-channel channel and the tap coef- 
30 ficient G(n) for linear filtering by the first round of adaptive equalization of said received signal R(n) of said 

M-channel, and outputting said first soft decision value X 1 [b(n)]: 
in the second and subsequent rounds of equalization, 

(b) calculating the likelihood b'(n) of a code bit sequence b(n) from said second soft decision value ^[b 
(n)] obtained by decoding said code bit b(n) based on said first soft decision value X..,[b(n)}; 

35 (c) calculating an estimated impulse response value vector H L (n) by approximated with the intersymbol 

interference component for the code bit b(n) regarded as zero, and calculating the tap coefficient G(n) for 
said linearfiltering from said estimated impulse response value vector H L (n) to update said tap coefficient 
G(n): 

(d) linear-filtering said likelihood b'(n) of said code bit sequence b(n) with said estimated impulse response 
40 value H L (n) to generate a replica; 

(e) subtracting said replica from said received signal R(n) to generate a difference signal R c (n) without 
intersymbol interference; 

(f) linear-filtering said difference signal R c (n) with said tap coefficient G(n) to generate a signal Z(n); and 

(g) outputting, as the results of said second and subsequent rounds of adaptive equalization, a first soft 
45 decision value X 1 [b(n)] updated with said signal 2(n) and said estimated impulse response value vector 

H L (n). 

3. The method of claim 2, wherein said step (a) comprises the steps of: 

50 (a-1) calculating the estimated impulse response value H m (n) of said M-channel channel based on said re- 

ceived signal R(n) and a known signal; 

(a-2) calculating said tap coefficient G(n) by an adaptive algorithm based on said received signal R(n) and 
said known signal; 

(a-3) linear-filtering said received signal R(n) with said tap coefficient G(n); and 
55 (a-4) calculating said first soft decision value A-,[b(n)] from said signal Z(n) and said estimated impulse response 

value H m (n). 

4. The method of claim 2 or 3, wherein said step (c) includes a step of calculating the tap coefficient G(n) for said 
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linear filtering from said estimated impulse response value vector H L (n) through utilization of a Matrix Inversion 
Lemma.. 



The method of claim 2, wherein said step (a) comprises the steps of: 



(a-1 ) calculating the estimated impulse response value H m (n) of said M-channel channel based on a sample 
value sequence of said received signal R(n) and a known signal; 

(a-2) determining from said estimated impulse response value H m (n) whether the received signal power is 
larger than a predetermined reference value, and, if larger than said reference value, storing the corresponding 
path as an effective path in a memory; 

(a-3) calculating said tap coefficient G(n) based on said received signal R(n) and said known signal, storing 
a tap coefficient G'(n) for said linear filtering corresponding to said effective path in a memory, and storing a 
received signal vector R*(n) corresponding to said effective path in a memory; and 

(a-4) calculating said first soft decision value A^fbfn)) from said estimated impulse response value H m (n), said 
tap coefficient G'(n) and said received signal vector R'(n). 



The method of claim 5, wherein: 



said step (c) comprises the steps of: 



(c-1) calculating the tap coefficient G'(n) for said linear filtering corresponding to said effective path from 
said estimated impulse response value H L *(n) composed of the component corresponding to said effective 
path through utilization of an Matrix Inversion Lemma; and 

(c-2) storing said tap coefficient G'(n) for said linear filtering corresponding to said effective path in a 
memory; 



said step (d) is a step of linear-filtering said tap coefficient G'(n) with said estimated impulse response value 
vector H L '(n) to obtain a replica signal; 

said step (e) is a step of storing in a memory a difference signal R c '(n) corresponding to said effective path, 
said difference signal R c '(n) being obtained by subtracting said replica signal from said received signal R'(n); 
said step (f) is a step of linear-filtering said difference signal R c *(n) with said rap coefficient G'(n) to generate 
a signal Z'(n); and 

said step (g) is a step of calculating said first soft decision value ^[bfn)] from said estimated impulse response 
value vector H L '(n) and said signal Z'(n). 

The method of claim 2, wherein, letting J represent the maximum number of delayed symbols to be considered, 
and letting a received signal sample vector of an M-channel channel be represented by r(n)=[r 0 (n)r 1 (n)...r M . 1 (n)] T , 
said received signal vector R(n) by R(n)=[r(n+J-1)r(n+J-2)...r(n)] T f a channel weighting coefficient vector by h(n; 
j)=[h 0 (n;j)h 1 (n;j)...h M . 1 (n;j)] T , and a channel matrix H m (n) of said estimated impulse response value by 



H m (n) = 



h(n;0) h(n;1) — h(n;J-l) 0 
0 h(n;0) h(n;l) b(n;J-l) 



0 
0 



0 
0 



h(n;0) h(n;l) 
0 h(n;0) 



h(n;J-J) 0 
b(n;l) — h(n;J-l) 



MJx(2J-l) 



said step (a) comprises the steps of: 



(a-1 ) linear-filtering a training signal b(n) in a training signal period with said estimated impulse response value 
H m (n) to generate a replica H m (n)b(n); 

(a-2) generating the difference between said received signal R(n) and said replica H m (n)b(n) as a difference 
vector R c (n); 

(a-3) linear-filtering said received signal R(n) with said tap coefficient G(n) to generate an output Z(n)=G(n) H R 

(n); 

(a-4) determining said tap coefficient G(n) by an adaptive algorithm based on the difference between said 
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10 



15 



20 



output Z(n) and said training signal b(n): and 

(a-5) calculating a soft decision value X 1 [b(n)]=4Real{Z(n)}/(1-u.) based on the estimated impulse response 
value H m (n) for linear filtering and said output Z(n), and outputting said soft decision value A^[b(n)]=4Real{Z 
(n)}/(1-p.) as the result of said first round of equalization; 

said step (b) is a step of calculating the likelihood b'(k)=tanh[X2[b(k)/2] of a code bit sequence b(k) from a soft 
decision value X^[b(n)] of a decoded bit provided as a priori information with said k set within the range of n-(j- 
1)<k<n+(J-1); 

said step (c) is a step of calculating said tap coefficient G(n) by approximating 

G(n)=[H L (n)H L (n) H (n)-a 2 ir 1 H L (n) 

H L (n)=[h 0 (n;J-1)...h^ 
said step (d) is a step of linear-filtering, with said H L (n), an estimated value vector, 

B , (nHb'(n+(J>1))b , (n+(J-2))...b , (n+1)Ob'(n-1)...b , (n-(J-1))3 T , 



of a code bit that affects, as intersymbol interference, said code bit b(n) at time n to thereby obtain a replica 
H L (n)B»(n); 

said step (e) is a step of calculating a difference vector R c (n)=R(n)-H L (n)B'(n) between said replica H L (n)B' 
(n) and said received signal R(n); 
25 said step (f) is a step of linear-filtering said difference vector R c (n) with said tap coefficient G(n) and outputting 

the result of said linear filtering Z(n)=G(n) H R c (n); and 
said step (g) is a step of obtaining a soft decision value 

30 Mb(n)J- ^(n) 

^(n)=H L (n) H G(n) 

35 as the output of said second and subsequent rounds of equalization from said output Z(n) and said estimated 

impulse response value vector H L (n). 

8. The method of claim 6, wherein, letting J represent the maximum number of delayed symbols to be considered, 
and letting a received signal sample vector of an M-channel channel be represented by r(n)=[r 0 (n)r 1 (n)...r M . 1 (n)) T , 
40 said received signal vector R(n) by R(n)=[r(n+J-1)r(n+J-2)...r(n)] T , a channel weighting coefficient vector by h(n; 

j)=[h 0 (n;j)h 1 (n;j)...h M _ 1 (n;j)] T , and a channel matrix H m (n) of said estimated impulse response value by 



45 



50 



H m (n) = 



fa(n;0) h(n;l) -« b(n;J-l) 0 
0 h(n;0) h(n;l) b<n;J-I) 



0 
0 



0 
0 



h(n;0) h(n;l) 
0 h(n;0) 



h(n;J-l) 0 
h(n;I) h(n;J-l) 



MJx(2J-l) 



55 



said step (a-3) comprises the steps of: 

(a-3-1) linear filtering a training signal b(n) in a training signal period with said estimated impulse response 
value H m (n) to generate a replica H m (n)b(n); 

(a-3-2) generating the difference between said received signal R(n) and said replica H m (n)b(n) as a difference 
vector R c (n); 
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10 



(a-3-3) linear filtering said received signal R(n) with said tap coefficient G(n) to generate an output Z(n)=G 
(n) H R(n); and 

(a-3-4) determining said tap coefficient G(n) by an adaptive algorithm based on the difference between said 
output 2(n) and said training signal b(n), and storing a received signal R'(n) and a tap coefficient G'(n) of those 
components of said received signal R(n) and said tap coefficient G(n) which correspond to said effective path : 
said step (a-4) is a step of calculating a soft decision value X 1 [b(n)]=4Real{Z(n)}/(1- u ,) based on the 
estimated impulse response value H m '(n) for said linear filtering and said output Z(n), and outputting said soft 
decision value A. 1 [b(n)]=4Real{Z(n)}/(1 -u.) as the result of said first round of equalization; 
said step (b) is a step of calculating the likelihood b'fkHanh^fbfkJte] of a code bit sequence b(k) from a soft 
decision value X^[b(n)] of a decoded bit provided as a priori information with said k set within the ranqe of n-(i- 
1)<k<n+(J-1); 

said step (c-1) is a step of calculating said tap coefficient G'(n) by approximating 
15 G , (nHH L '(n)H L '(n) H (n)-o 2 l]" 1 H L , (n) 

H L '(nMh 0 (n;J-1)...h^ 

20 sa 'd step (d) is a step of linear filtering, with said HL'(n), an estimated value vector 

B , (nHb , (n+(J-1))b»(n+(J-2))...b'(n+1)0b'(n-1)...b , (n-{J-1))] T > 



25 



30 



35 



of a code bit that affects, as intersymbol interference, said code bit b(n) at time n to thereby obtain a replica 
H L (n)B'(n); 

said step (e) is a step of calculating a difference vector R c .(n)=R'(n)-H L (n)B'(n) between said replica H L '(n)B' 
(n) and said received signal R'(n); 

said step (f) is a step of linear-filtering said difference vector R c '(n) with said tap coefficient G'(n) and outputting 
the result of said linear filtering Z'(n)=G , (n) H R c '(n); and 
said step (g) is a step of obtaining a soft decision value 



^^ 4 T y £ <n)} 

^(n)=H L '(n) H G'(n) 

as the output of said second and subsequent rounds of equalization from said output Z'(n) and said estimated 
40 impulse response value vector H L '(n). 

9. An adaptive equalizer comprising: 

an impulse response estimating part for calculating an estimated impulse response value of each of Channel 
45 based on a received signal and a known signal; 

a tap coefficient calculating part for calculating the tap coefficient of a linear filter by an adaptive algorithm 
based on said received signal and said known signal; 

said linear filter having set therein said tap coefficient, for linear-filtering said received signal; and 
a soft decision value calculating part for calculating a soft decision value from said estimated impulse response 
50 value and the result of sard linear filtering. 

10. The equalizer of claim 9 further comprising: 

a storage part for storing said estimated impulse response value; 
55 a likelihood calculating part for calculating, from said soft decision value, its likelihood; and 

means for obtaining a replica signal by subjecting said likelihood to linear filtering with an estimated impulse 
response value vector obtained from said stored estimated impulse response value vector by approximating 
intersymbol interference with the code b(n) concerned to zero; 
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10 



wherein said tap coefficient calculating part includes means for calculating the tap coefficient of said linear 
filter by an adaptive algorithm based on said received signal and said known signal when a soft decision value is 
input from a decoder to said tap coefficient calculating part, and for calculating said tap coefficient from said esti- 
mated impulse response value vector when said soft decision value is input to said tap coefficient calculating part; 
and 

said linear filter performs linear filtering of said received signal in the absence of said soft decision value and 
the difference between said received signal and said replica signal through the use of said tap coefficient in the 
presence of said soft decision value, and provides the linear filtering output to said soft decision value calculating 
part. 

11. The adaptive equalizer of claim 9, further comprising: 



a path decision part for determining from said estimated impulse response value of said each channel whether 
the received power of the corresponding path is larger than a predetermined reference value;; 
is a path memory for storing as an effective path a path determined as being larger than said predetermined 

reference value; 

a tap coefficient memory for storing as a new tap coefficient the component of that one of said tap coefficients 
corresponding to said effective path; and 

a signal identifying part for identifying that received signal component of said received signal corresponding 
20 to said effective path: 

wherein said soft decision value calculating part calculates said soft decision value from the received signal 
component corresponding to said effective path : the estimated impulse response value corresponding to said 
effective path and the tap coefficient corresponding to said effective path. 



25 



12. The adaptive equalizer of claim 11 s further comprising: 



a likelihood calculating part for calculating the likelihood of a code from said soft decision value; 
a replica generating linear filter for generating a replica signal of said received signal by linear-filtering said 
30 likelihood with an estimated impulse response value composed of that component of said estimated response 

value vector corresponding to said effective path; a subtractor for subtracting said replica signal from that 
component of said received signal corresponding to said effective path to obtain a difference signal; and 
a difference memory for storing only that difference signal component of said difference signal corresponding 
to said effective path: 



35 



wherein sad tap coefficient calculating part calculates, in second and subsequent rounds of equalization, the 
tap coefficient of said linear filter from said estimated impulse response value vector corresponding to said effective 
path through the use of a Matrix Inversion Lemma: and 

said soft decision calculating part is a means for calculating, in said second and subsequent rounds of equal- 
ization, said soft decision value from said estimated impulse response value vector corresponding to said effective 
path, said tap coefficient of said linear filter corresponding to said effective path and said difference signal corre- 
sponding to said effective path. 



45 



50 



55 



20 



WSDOCIO: <EP 121 7798 A2_l_> 



EP 1 217 798 A2 




21 



BNSDOCID: <EP 1217798A2_I_> 



EP 1 217 798 A2 



J- 



31 



JL 



31 



31 

















1 I 




r " n 




< 



COEFFICIENTS 

ADAPTATION 

CIRCUIT 



34 



7 



OUTPUT 



32 

VARIABLE 

WEIGHTING 

CIRCUIT 



FIG. 2 



3NSDOCID: <EP 1217798A2_I_> 



22 



EP1 217 798 A2 



112 



113- 



ix £ 



~y 110 LINEAR FILTER PART 



TAP 
COEFFICIENT 
CALCU 



160 

X 



G(n) n 



111 



Z(n) 



TRAN 
FILTER 



1 



120 



SOFT DECISION 
VALUE CALCU 



j.j 



104 



R(n) 



TRAN 
FILTER 



■fc 



105 



ik H m (n) 



•103 



b(n) 



/ 
100 

IMPULSE RESPONSE 
ESTIMATING PART 



CALCU 
PART 


^-106A 


; STORAGE 
• PART 


^106B 


. t 




Rc(n) 


-I 

106 



Ai(b(n)] 



101 

> 



TRAINING 
SIGNAL 
102 



V 

130 
LIKELIHOOD 
CALCU 



150 

DISCRIMINATING 
PART 



ADAP ESTIMATING 
PART 



FIG. 3 



BNSDOCID: <EP 1217798A2_I_> 



23 



EP1 217 798 A2 



112 



113 



H L (n) 



TAP 

COEFFICIENT 
CALCU 



f G(nU 



, ^ 110 LINEAR FILTER PART 

v-f 



111 



L 



160 

\ 

Q 



Z(n)j 



104 

I. 



R(n) 





TRAN 




FILTER 








TRAN 




FILTER " 
■ 177 



120 



SOFT DECISION 
VALUE CALCU 



%(n) 



103 



H m (n) 



|b'(n) 



Ai[b(n)} 
► 



101 



[ CALCU 1 


^106A 


i PART 


r 
t 


STORAGE 


^-106B 


PART 



105 _ [[PART K 1 



TRAINING 
SIGNAL 
^102 

131 



X 2 [b(n)] 



T 
130 
LIKELIHOOD 
CALCU 

, 100 
I IMPULSE RESPONSE 
! ESTIMATING PART 



i 

k 



106 

ADAP ESTIMATING 
PART 



FIG. 4 



INSDOCID: <EP 1217798A2_I_> 



24 



EP 1 217 798 A2 



LU 
CD 

LU 

O 

LU 

> 



10- 1 

10- 2 
103 

10- 5 

10- 6 





















- 










































\% 
\* 
\ % 
\ * 
\ % 

\ V 

\ » 


















3 




2 3 4 5 6 
AVERAGE Eb / No (dB) 

POWER 




d 1 d 2 



— o— AO 
-b-A1 
- BO 

-13- B1 



8 



FIG. 5 



d 3 d 4 DELAY TIME 



FIG. 6 



25 



BNSDOCID: <EP 121 7798A2_I_> 



EP 1 217 798 A2 



TAP 

112 COEFFICIENT 



"~l 



CALCU 



113 



r 



110 LINEAR FILTER PART 



CALCU 
PART 


^112A 
^112C 


ID PART 


STORAGE 
PART 


^112B 



111 



TRAN 
FILTER 



TRAN 
FILTER 




+ • 



Z-(n) 



Z(n) 



.120 



SOFT DECISION 
VALUE CALCU 



-103 



ID PART 
140 



H m (n) 



HJn) 



1 



!b(n) 



101 



CALCU 
PART 


^-106A 
^106C 

^106B 


ID PART 


STORAGE 
PART 





TRAINING 
SIGNAL 
102 

-6 



R c (n) 



V 
130 
LIKELIHOOD 
CALCU 



I 100 

j IMPULSE RESPONSE 
. I ESTIMATING PART 



106 

ADAP ESTIMATING 
PART 



FIG. 7 



BNSDOCID: <EP 1217798A2J_> 



26 



EP 1 217 798 A2 



r TAP 

H2 COEFFICIENT 113 
nCALCU 



£ 



CALCU t<112A 
PART v 

L<112C 



ID PART 



STORAGE ^H2B 
PART * 



L. 



160 



\ 

Q 



111 



+ ! 



■ 

r 



110 LINEAR FILTER PART 



Z(n) 



Z(n): 



1 



.120 



SOFT DECISION 
VALUE CALCU 



TRAN 
FILTER 



ID PART 
140-^ 



170. 



104 

/ 
O 



R(n) 



TRAN 
LTER 



105 



103 



H L (n) 



H L '(n) 
b'(n) 



CALCU 
PART 



ID PART 



STORAGE 
PART 



106A 
106C 

•106B 



7 



Aitb(n)] 
► 



101 



TRAINING 
SIGNAL 
102 

4 



T 

130 
LIKELIHOOD 
CALCU 



Rc'(n) [ J 



100 

IMPULSE RESPONSE 
ESTIMATING PART 



106 

ADAP ESTIMATING 
PART 



FIG. 8 



27 

BNSDOCID: <EP 1217798A2_I_> 



EP 1 217 798 A2 



POWER 



EQUALIZATION RANGE 



(MAXIMUM NUMBER OF DELAYED 
SYMBOLS TO BE CONSIDERED IN 
ADAPTIVE EQUALIZER) 



Liill 



lllllll 



DELAY TIME 

FIG. 9 



POWER 



THRESHOLD 
VALUE 



EQUALIZATION RANGE 



(MAXIMUM NUMBER OF DELAYED 
SYMBOLS TO BE CONSIDERED IN 
q ADAPTIVE EQUALIZER) 

° 0 O: EFFECTIVE PATH 

Q TQ 



Mi l l 



lllllll 



DELAY TIME 



FIG. 10 



28 



NSDOCID: <EP 121 7798 A2_l_> 



EP 1 217 798 A2 




POWER 



0 1T 2T 3T 4T DELAY TIME 



FIG. 12 



29 



BNSDOCID: <EP 1217798A2J_> 



J 




Europaisches Patentamt 
(19) ^j))) European Patent Office 

Office europeen des brevets EP 1 21 7 798 A3 

(12) EUROPEAN PATENT APPLICATION 

(88) Date of publication A3: (51) , nt c , ?, H0 4L 25/03, H04L 25/02 

14.04.2004 Bui.etin 2004/16 H04L 1/00, H03M 13/29 ' 

(43) Date of publication A2: 

26.06.2002 Bulletin 2002/26 



(21) Application number: 01129388.3 

(22) Date of filing: 18.12.2001 



(84) Designated Contracting States: 

AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 
MC NL PT SE TR 
Designated Extension States: 
AL LT LV MK RO SI 

(30) Priority: 19.12.2000 JP 2000385699 

(71) Applicant: NTT DoCoMo, Inc. 
Tokyo 100-6150 (JP) 

(72) Inventors: 

• Ohmori, Hiroo 

Yokohama-shi, Kanagawa 236-0031 (JP) 



(54) 



• Asai, Takahiro 

Yokosuka-shi, Kanagawa 239-0806 (JP) 

• Matsumoto, Tadashi 

Ebina-shi, Kanagawa 243-0431 (JP) 

• Tomisato, Shlgeru 
Kanagawa 233-0008 (JP) 

(74) Representative: Hoffmann, Eckart, Dipl.-lng. 
Patentanwalt, 
Bahnhofstrasse 103 
82166 Grafelflng (DE) 



Adaptive equalization method and adaptive equalizer 



CO 
< 
00 



(57) During reception of a training signal b(n) in a 
received signal R(n) an estimated impulse response H m 
(n) of an M-channel channel and a tap coefficient G(n) 
of a linear filter (111) are calculated by an adaptive al- 
gorithm through the use of the received signal R(n) and 
the training signal b(n). For an information symbol of the 
received signal R(n), the received signal R(n) is subject- 
ed to linear filtering with the most recently calculated tap 
coefficients G(n), and the linear filtering output Z(n) and 
the most recently estimated impulse response value H m 
(n) are used to calculate a soft decision value A. v In the 
second and subsequent rounds of equalization., the like- 
lihood b'(n) of a soft decision value A^bfn)] from a de- 
coder is calculated, and a replica are generated by lin- 
early filtering the likelihood b'(n) with an estimated im- 
pulse response value vector H L (n) obtained by approx- 
imating intersymbol interference with the current code 
b(n) to zero. A difference signal R c (n) between the rep- 
lica and the received signal is calculated, and the esti- 
mated impulse response value vector H L (n) is used to 
update the tap coefficients G(n). Then the signal Z(n) is 
obtained by linearly filtering the difference signal R c (n) 
with the updated tap coefficient G(n), and the signal Z 
(n) and the estimated impulse response value vector H L 
(n) are used to calculate the soft decision value a. v 



~Y 1 10 LINEAR FILTER PART 




/ 
100 

IMPULSE RESPONSE 
ESTIMATING PART 



106 

ADAP ESTIMATING 
PART 



FIG. 3 



CL 
LU 



Printed by Jouve, 75001 PARIS (FR) 



3NSDOCID: <EP 1217798A3_I_> 



EP 1 217 798 A3 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 01 12 9388 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (lnt.CI.7) 



D,P, 
X 



NIKOLAI NEFEDOV; MARKKU PUKK1LA: "Turbo 
Equalization and Iterative (Turbo) 
Estimation Techniques for Packet Data 
Transmission 11 

2ND INTERNATIONAL SYMPOSIUM ON TURBO CODES 

AND RELATED TOPICS, 

4 September 2080 (2080-09-04), - 7 

September 2000 (2000-09-07) pages 423-426, 

XP0O09O25553 

* page 424, paragraph 3.1 - page 425, 
paragraph 4 * 

* figure 1 * 

WANG X ET AL: "INTERACTIVE (TURBO) SOFT 
INTERFERENCE CANCELLATION AND DECODING FOR 
CODED CDMA" 

IEEE TRANSACTIONS ON COMMUNICATIONS, IEEE 
INC. NEW YORK, US, 

vol. 47, no. 7, July 1999 (1999-07), pages 
1046-1061, XP000849079 
ISSN: 0090-6778 

* page 1047, paragraph II -A - page 1048, 
paragraph III * 

* page 1050, paragraph IV. B 

* page 1055, paragraph V.A 
paragraph V.B * 

* page 1057, paragraph V.C 
paragraph V.D * 



1.9 



H04L25/03 
H04L25/02 
H04L1/00 
H03M13/29 



1-12 



TECHNICAL FIELDS 
SEARCHED (lnLCI.7) 



- page 1052 * 
page 1056, 

page 1058, 



H03M 
H04L 



REYNOLDS D ET AL: n Low -complexity 

Turbo-equalization for diversity channels" 

SIGNAL PROCESSING, AMSTERDAM, NL, 

vol. 81, no. 5, May 2001 (2001-05), pages 

989-995, XP004239259 

ISSN: 0165-1684 

* the whole document * 



1-12 



-/- 



The present search report has been drawn up for ail claims 



Place o< search 

Munich 



Date of completion of the search 

12 February 2004 



Examiner 

Schiffer, A 



CATEGORY OF CITED DOCUMENTS 

X : particuia/ty relevant if taken alone 

Y ; particularly relevant ri combined with another 

document ot the same category 
A : technological background 
O : non-written disclosure 
P : intermediate document 



T : theory or principle underlying the invention 
E : earlier patent document but published on. or 

after the f if ng dale 
D : document cited in the application 
L : document cited for other reasons 



& : member of the tame patent family, corresponding 



2 



3NSDOCID: <EP 1217798A3J_> 



EP 1 217 798 A3 




European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 01 12 9388 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (lnt.CL7) 



SELLATHURAI M ET AL: "TURBO-BLAST for 
high-speed wireless comnunications" 
PROCEEDINGS OF IEEE CONFERENCE ON WIRELESS 
COMMUNICATIONS AND NETWORKING, 
vol. 1, 23 September 2000 (2000-09-23), 
pages 315-320, XP010532516 
PISCATAWAY, NJ, IEEE, USA 
* figure 2 * 



U9 



TECHNICAL FIELDS 
SEARCHED (lnt.CL7) 



The present search report has been drawn up for all claims 



Plac* of teat ch 


DaJ« c4 completion ol the search 


Examiner 


Munich 


12 February 2004 


Schiffer, A 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant t taken alone 

Y : particularly relevant if combined wfth another 

document of the tame category 
A : technological background 
O : non written disclosure 
P . intermediate document 



T : theory or principle underlying the invention 
E : earftet patent document, but pubfcshed on, or 

after the fi&ng date 
O : document cited in the application 
L : document cited for other reasons 



& : member of the some patent family, corresponding 
document 



3 



BNSDOCID: <EP 1217798A3_I_> 



5 



THIS PA61 BLANK <usero) 



